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The Nariai black hole, whose two horizons are lying close to each other, is an extreme and 
important case in the research of black hole. In this paper we study the evolution of a massless 
scalar field scattered around in 5D Schwarzschild-de Sitter black string space. Using the method 
shown by Brevik and Simonsen (2001) we solve the scalar field equation as a boundary value problem, 
where real boundary condition is employed. Then with convenient replacement of the 5D continuous 
potential by square barrier, the reflection and transmission coefficients (R, T) are obtained. At last, 
we also compare the coefficients with usual 4D counterpart. 
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I. INTRODUCTION 



Black hole radiation, which was proved originally by Stephen Hawking with the method of quantum field on 
gravitational collapsing [l[ 0, indicates that black holes are not perfect black but radiate thermally and eventually 
explode. Since then, many people have used various methods and techniques to research black hole through the 
^' j ' particles radiating from it, such as the simple Klein-Gordon particles and Dirac particles (for some early works, see 
Damour and Ruffini @] and Chandrasekhar Q respectively). Here, scalar particles are only considered. Recently, 
Higuchi et al. Q and Grispino et al. J6| gave the scalar field solution outside a Schwarzschild black hole; Brady et al. 
| Q, Brevik et al. Q and Tian et al. [!| studied the Schwarzschild-de Sitter (SdS) case; Guo et al. [lCj made further 
studies in the Reissner-Nordstrom-de Sitter one. As for recent studies of searching evaporating black holes, one can 
O refer to the works [H- 

■ The idea that the world may have more than four dimensions is due to Kaluza [l4| and Klein [TH, who realized 
that a 5D manifold could be used to unify general relativity with Maxwell's theory of electromagnetism. After that, 
many people focus on the robust higher dimensional space. Here, we consider the Space- Time-Matter (STM) theory 
t**"*- 1 presented by Wesson and co-workers [l|| ^7|- This theory is distinguished from the classical Kaluza-Klein theory 
for a non-compact fifth dimension, the 4D source is induced from an empty 5D manifold. Because of this, the STM 
theory is also called induced matter theory and the effective 4D matter is called induced matter. That is, in STM 
theory, 5D manifold is Ricci-flat while 4D hypersurface is curved by the 4D induced matter. Mathematically, this 
approach is supported by Campbell's theorem which states that any analytical solution of N-dimensional Einstein 
equations with a source can be locally embedded in an (N-fl)-dimensional Ricci-flat manifold [l8|. In the framework 
of STM, people studied many works such as Quantum Dirac Equation [l9[ , Perihelion Problem [2(| , Kaluza-Klein 
Solitons [2l|, Black Hole [H J24J, Solar System Tests [H| and so on. 

In order to avoid interactions beyond any acceptable phenomenological limits, people assume standard model fields 
(such as fermions, gauge bosons, Higgs fields) are confined on a (3 + 1) dimensional hypersurface (3-brane) without 
accessing along the transverse dimensions. The branes are embedded in the higher dimensional spacetime (bulk), in 
which only gravitons and scalar particles without charges could propagate under standard model gauge group. There 
are also man y w orks (for a review with large extra dimensions see [26]) focusing on Hawking radiation such as (27j 

[H H [H in. 

Cosmological constant A, which is a parameter with dimension L~ 2 (L is length), is one of focuses in Gravitation 
Theory. The acceleration of the universe is explained by the required repulsive force produced by a non-zero and 
positive cosmolo gica l constant A. Current Snla observation data shows that the cosmological constant has a value of 
A ~ 10~ 52 m~ 2 [32| (33| [HI- Its robust non-zero magnitude engages the researching interest in the space contained 
cosmological constant. Especially the black hole contained effective cosmological constant is studied widely either in 
higher dimensions background [2|| (2?| [Hj] or in usual 4D case Q [U @ [HI • Sometimes for the sake of study, A 
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is considered as a free parameter like in the works [8| [9( [lOj]. In SdS space, the interval between black hole horizon 
r e and cosmological horizon r c becomes smaller with increase value of A. If cosmological constant A reaches its 
maximum, Nariai black hole will be arisen. In this paper, we study how a massless scalar field evolves in this extreme 
case. 

This paper is organized as follows: in section II, the 5D SdS black string space, the time-dependent radial equation 
about Rcjir, t) and the fifth dimensional equation about L(y) are restated. In section III, by a tortoise coordinate 
transformation, the radial equation becomes a Schrodinger wavelike one. According to the boundary condition and the 
tangent approximation, a full numerical solution is presented. In section VI, using the replacement of real potential 
barriers around black hole by square barriers, the reflection and transmission coefficients are obtained. Section V is 
a conclusion. 

We adopt the signature (+, -, -, -, -), put K, c, and G equal to unity. Lowercase Greek indices /i, v, . . . will be taken 
to run over 0, 1, 2, 3 as usual, while capital indices A, B, C . . . run over all five coordinates (0,1,2,3,4). 



II. THE MASSLESS SCALAR FIELD IN 5D SCHWARZSCHILD-DE SITTER BLACK STRING SPACE 

Within the framework of STM theory, an exact 5D solution presented by Mashhoon, Wesson and Liu [l(| [Hf [36| 
describes a 5D black hole. The line element takes the form 



dS 2 = 



Ae 



f(r)dt 2 



1 



■dr 2 - 



sin' 



: 9d(/) 2 ) 



-de 



In our case 



tl N 1 2M A 2 

f(r) = 1 - — - gr 2 , 



(1) 



(2) 



where £ is the open non-compact extra dimension coordinate, A is the induced cosmological constant and M is the 
central mass. The part of this metric inside the square bracket is exactly the same line-element as the 4D SdS solution, 
which is bounded by two horizons — an inner horizon (black hole horizon) and an outer horizon (one may call it 
cosmological horizon). This metric (pj satisfies the 5D vacuum equation Rab = 0, therefore, there is no cosmological 
constant when viewed from 5D. However when viewed from 4D, there is an effective cosmological constant A. So one 
can actually treat this A as a parameter which comes from the fifth dimension. This solution has been studied in 
many works [37| focusing mainly on the induced constant A, the extra force and so on. 
We redefine the fifth dimension in this model, 



(3) 



Then we use dTJ ~ (J3J) to build up a RS type brane model in which one brane is at y = 0, and the other brane 
is at y = Hence the fifth dimension becomes finite. It could be very small as RS I brane model [38| or very 
large as RS II model [3!|. The relation between STM theory and brane world theories, and the embedding of 5D 
solutions to brane models are studies in [4(| [42j [Hj]. For the present brane model, when viewed from a (£ or 
y = constant) hypersurface, the 4D line-element represents exactly the SdS black hole. However, when viewed from 
5D, the horizon does not form a 4D sphere — it looks like a black string lying along the fifth dimension. Usually, 
people call the solution to the 5D equation ^Gab — A^qab (As is the 5D cosmological constant) as the 5D SdS 
solution. Therefore, to distinguish with it, we call the solution ([T]) a black string, or more precisely, a 5D Ricci-flat 
SdS solution. 

After redefining the fifth dimension, the metric |T]) can be rewritten as 



dS 2 = eV¥» 



f{r)dt 2 



dr 2 



(de 2 



'■) - dy< 



where y is the new fifth dimension. Expression Q can be recomposed as follows 

f(r) = — (r - r e )(r c - r)(r - r Q ). 
3r 



(4) 



(5) 



The singularity of the metric ((H) is determined by f(r) — 0. Here we only consider the real solutions. The solutions 
to this equation are black hole event horizon r e , cosmological horizon r c and a negative solution r = — (r e + r c ). The 
last one has no physical significance, and r c and r e are given as 



3 



?7- 



2 

(6) 



rc = ^cosr;, 
= ^ cos(120° - 77) 



where ^ = 5 arccos(— 3MVA) with 30° < rj < 60°. The real physical solutions are accepted only if A satisfy AM 2 < i 

& 

Then we consider a massless scalar field $ in the 5D black string spacetime, obeying the Klein-Gordon equation 

□$ = 0, (7) 

where □ = -^g§w {y/99 AB ~b§b) ^ s the 5D d'Alembertian operator. We suppose that the separable solutions to Eq. 
([7]) are in the form 

* = -jL=-R^{r, t)L(y)Y lm (6, 0), (8) 

where Ruj{r, t) is the radial time-dependent function, Yi m (Q,<p) is the usual spherical harmonic function, and L(y) is 
the function about the fifth dimension. The differential equations about y and (t, r) are 

d 2 L(y) , FKdUy) nr , , „ 



dy 2 V 3 dy 

-r 



fk'i (t) + S HI (t)) - ^ + i « + ^ T " °- < 10 > 

Eq. (|10[) is a time-dependent radial differential equation. Eq. ^ is a differential equation about y, where Q is a 
constant which is adopted to separate variables (t, r, 0, 0, y). 



III. THE NARIAI BLACK HOLE AND ITS BOUNDARY VALUE PROBLEM 

Nariai black hole [44| [45[ occurs when the cosmological horizon is very close to the black hole horizon r c — ► r e . It 
is an extreme and important kind of SdS black holes. The cosmological constant in this limit is given by 

KM 2 = -. (11) 
9 V 1 

Substituting Eq. (fTTjl into Eq. ((6]), we can get 77 = 60° and the horizons = r e = r c = 3M. As an illustration of 
the accuracy , we mention that the choice AM 2 = 0.11 Q @ leads to r e = 2.8391M and r c — 3.1878A/. In order to 
simplify numerical calculation, we will put M = 1 in this paper. 



A. The Fifth Dimensional Function L(y) 

In our previous paper [24j], we have introduced a massless scalar field to stabilizing this black string brane model. 
Considering a single mode of the scalar field, the wave function for this mode may reach its maximum value but 
keep smooth and finite at the brane. Hence, a steady standing wave is constructed. A suitable superposition of 
some of the quantized and continuous components of L(y) may provide a wave function which is very large at y = 
and drops rapidly for y ^ 0. Naturally, a practical 3-brane is formed at the y — hypersurface. According to this 
"standing wave" condition in the bulk, the spectrum of Q is broken into two parts. One is the continuous spectra 
below I A and the other is the discrete spectra above |A. The quantum parameter f2„ is 

3 

where n=l, 2,3... and yi is the thickness of the bulk. So the solutions to Eq. ([9]) are 



n n = + -a, (12) 



Ce-^y cos (nn^j , n = 1,2,3- ••, for Q > 
L (V) = < (Ci + C 2 y) e-^, for n = M (13) 



3A 
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FIG. 1: The first three eigenfunctions of L n (y) in Nariai case: Li(solid), Z/2(dashed), and /^(dotted) with 3/1 — > 1, C=l and 



where A = 0.11 is the cosmological constant, and (C, yo), (Ci, C2), (C3, C4) are the three pairs of integration constants. 
As an illuminating example, Fig. [T]depicts the quantized sates of L(y). It illustrate that the eigenfunctions get the 

maximum on the brane y = and get an extremum on the other brane y = y\. For an exponential factor e~^^ y 
in the first solution of Eqs. (fT3|) . the wave function L n decays along the fifth dimension. Comparing with general 
result [24}, extreme cosmological constant A = 0.11 gives a more fiercely decay solutions. As the fifth dimension 
becomes bigger, the probability |L(y)| 2 deflects from the original value (\L(y)\y =0 = 1) more and more larger. It gets 
an extremum instead of the maximum value on the other brane. In this way, two branes can be stabilized by scalar 



A more important aspect of scalar field is radial direction. In Eq. (|10p time variable can be eliminated by the 
Fourier component e~' luJt via 



where the subscript n presents a new wave function unlike the usual 4D case tp u i . So Eq. (|10p can be rewritten as 



A = 0.11. 



field. 



B. The Schrodinger wavelike equation 



(14) 



(15) 



where the potential function is given by 



V(r) 



1 df(r) 1(1 + 1) 



(16) 



r dr r 2 



Now we introduce the tortoise coordinate 




(17) 



The tortoise coordinate can be expressed by the gravitation surface as follows 
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where 



Explicitly, we have 
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A, 
A, 
■A. 



So under the tortoise coordinate transformation (|17p . the radial equation (|15[) can be rewritten as 

dx 2 



■4MV(r) 



Vuln(x) = 4M 2 w 2 $ u! „(a;) 1 



(19) 

(20) 
(21) 
(22) 

(23) 



which likes the form of Schrodinger equation in quantum mechanics. Notice that there are two various coordinates — 
r and x in this equation. So people also call it Schrodinger wavelike equation. The incoming or outgoing particle flow 
between event horizon r e and cosmological horizon r c is reflected and transmitted by the potential V(r). Substituting 
quantum parameters f2„ (jl2[) into Eq. (|16[) . the quantum potentials are obtained as follows 



V n (r) = f(r) 



1#) | 1(1 + 1) | r»V 
r dr r 2 y\ 
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(24) 



It is highly localizing near r ~ (r e + r c )/2 ~ 3, falling off exponentially in x at both r = r e and r — r c . Comparing 



with the similar case of usual 4D 
form of the potential for n=l, 2, 3 are illustrated in Fig. [21 



two additional monomials, 



and ^A, have appeared in the potential. The 




FIG. 2: The potentials of Nariai black hole with n=l (solid), n=2 (dotted), and n=3 (dashed). Here M=l, A = 0.11, / = 1, 
yi = 10 3//2 (a very large 5th dimension). The usual 4D potential (fi = 0) is also plotted with dash-dot line for comparison. 
The black hole horizon locates at the point Ar e ~ 2.8391 and the cosmological horizon locates at the point B. The potential 
tends to zero exponentially quickly as x — » ±oo. 



The numerical solution 



Near the horizons r e and r c , x — ► ±oo. According to Eq. ([5]) and Eq. we can get 

V(r e ) = V(r c ) = 0. (25) 



-1 



-40 -20 20 40 

x 

FIG. 3: Variation of the field amplitude versus x with M=l, 1=1, A = 0.11, yi = 10 3,/2 and n=l. The solution is close to a 
harmonic wave. 



So Eq. 



reduces to 



dx 2 



AM 2 u 2 



Vu,ln(x) = 0. 



Its solutions are e ± l2Mwx . i n this paper, we only take into account real field and choose the solution 

^uin = cos(2Mwie) 



(26) 



(27) 



as boundary condition near the two horizons (|6j). 

In real scalar field case, there are two methods to solve Schrodinger wavelike equation (|23p . One is tangent 
approximation 0] and the other is polynomial approximation 0. With any assigned cosmological constant A, we can 
always find an appropriate approximate method from those only by adjusting the parameters. Here the former one is 
adopted to analyze this model. Hence, we employe AM 2 = 0.11 and use the useful tangent approximation Q 



x(r) = 15tan[6(r - d) + 5], 

in which b = 2.7/ (r c — r e ) and d — (r c + r e )/2. Because the approximation ([28]) does not allow 
large, we shorten the interval of x to [-100,100]. So boundary condition (|27|) is rewritten as 

*»in(-100) = ^wlnClOO) = cos(200Mw). 



(28) 

to become very 



(29) 

Considering boundary condition (f2"9"]) and tangent approximation (|2"5)) . we can solve Eq. (f2"3")) numerically as a boundary 
value problem by Mathematica software. The variation amplitude of waves ty^i versus tortoise coordinate x is 
illustrated in Fig. [3] Considering actual circumstance, we use tortoise transformation (|18[) and also plot the amplitude 
versus r in Fig. [51 where we only give the first quantum state (n=l). The others can be treated by the same way. 



IV. THE REFLECTION AND TRANSMISSION 



We assume that the particle flux with energy E bursts towards a square well along the positive direction of x axis, 
where the potential is 

x ' \ 0, x < x\ or x > X2- 

From the view of quantum mechanics, considering the wave behavior of the particles, this process is similar to 
scattering on the surface of propagation medium with thickness of \x% — xi\. Parts of them are transmitted and parts 
of them are reflected back. According to statistical interpretation of wave function, whether the energy E > Vo or 
not, there is definite probabilities to transmit or reflect by the potential. The reflection and transmission coefficients 
denote the magnitude of those probabilities. 
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FIG. 4: Variation of the field amplitude versus r with M=l, 1=1, A = 0.11, yi = 10 3/2 and n=l. The waves pile up near the 
horizons. 



As mentioned above, it is necessary to replace the continuously varying potential barrier with a discontinuous barrier 
of constant height in analytical work. Therefore, the usual reflection and transmission coefficients can be obtained. 
With the method of [48| and Q, We suppose a scalar wave propagates from — oo to +oo, which is illustrated in Fig. 
O The same denotation is cited here, namely associating "1" with the incoming wave in the region — oo < x < x\, 
"2" with the potential plateau X\ < x < x%, and "3" with the outgoing wave in the region x 2 < x < +oo. Hence, 
potential V(x) in Eq. (f23|) reduces to 




According to square barrier §ST 



V(x) 

the solutions to Eq 



-oo < X < Xi, 
Xi < X < X2, 
X2 < X < +O0. 



(31) 



are 



ikxx _|_ fo^-ikxx 

^Loin — \ a 2 e + o 2 e , 



a 3 e 



— OO < X < Xl, 
Xl < X < x 2 , 
X 2 < X < +O0, 



(32) 



where ki = y 4M 2 (w 2 — Vi) (i — 1,2,3) are the wave numbers; ai and hi are the undetermined coefficients to the 
solutions . Then we define reflection coefficients for the plane interfaces dividing two media 



Rij — 



1 - Zi, 



1 



(33) 



where Zi, 



yy — -A- are the real impedance ratios between medium i and j. The width of barrier is d — x 2 — x\ and the 

height of square barrier is H = V 2 . So in this model reflection coefficients R and transmission coefficients T are given 
as 



R 



T = 



R12 + R23 + 2^/R 12 R 23 cos(2fc 2 cQ 
1 + Ri 2 R 23 + 2^R 12 R 23 cos(2fc 2 d) ' 
1 



1G 



(1 + Z 12 ) 2 {1 + Z 23 ) 2 1 + R 12 R 23 + 2^R 12 R 23 cos(2fc 2 d) ' 



(34) 
(35) 



Because the same width are adopted here, we only give the functional image of logi? versus hight H (or V 2 ) in Fig. 
[5l There is no surprise that it take oscillating like forms. One can read this feature directly from Eq. (f3~4"|) , which 
contains cosine functions. Then we use tangent approximation (|28[) and get the replacements of the 5D continuous 
potentials by square barriers in Fig. [6l The different reasonable V 2 are read off as the height of those square barriers. 
Meanwhile, we choose the incoming wave number to be ki — 2 (V\ — 0) and k 3 = 2 (V 3 = 0). Substituting those 
parameters into Eqs. (|34|) (|35[) . we can obtain reflection and transmission coefficients. Comparing those coefficients 
with usual 4D SdS case [8|, one can see the difference clearly in Table Q] Viewing from Fig. [51 we can get a relationship 
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FIG. 5: log R versus height H (or V 2 ) with the width d = 40, M = 1, I = 1, A = 0.11 and j/i = 10 3/2 
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FIG. 6: Replacement of real 5D SdS potential barriers around Nariai black hole by square barriers with n=l (solid), n—2 
(dotted), n=3 (dashed) and £7 = (usual 4D case with dash-dot line). We use M=l, 1=1, A = 0.11 and Vl = 10 3/2 . 

of the four heights Hq =0 < H n =\ < H n=2 < H n = 3 < H (the horizontal ordinate of the first extreme point). Hence, 
we can say that the four modes (£1 = 0, n = 1, n = 2, n = 3) are in the same monotone increasing space. Obviously, 
the reflection coefficients R (or T, notice R + T = 1) of 5D SdS black string are bigger (or smaller) than 4D case. 
Else, i2| n= i < R\ n=2 < R\ n =3- 



TABLE I: The reflection and transmission coefficients 



mode 


Xl 


%2 


d 


vi (or H) 


R 




4D SdS 


-20 


20 


40 


7.3 x 10~ 4 


1.3 x 10" 
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n=l 


-20 


20 


40 


10 x 10~ 4 


2.5 x 10" 
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n=2 


-20 


20 


40 


11 x 10~ 4 


3.0 x 10" 
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n=3 


-20 


20 


40 


12.5 x 10~ 4 


3.9 x 10" 
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V. CONCLUSION 

In this paper we have solved the real scalar field $ and obtained reflection and transmission coefficients (R, T) 
around the Nariai black hole in the 5D SdS black string space. We summarize what have been achieved. 

1. The 5D solution presented by Mashhoon, Wesson and Liu [3(| (3f| [l(| is exact in higher dimensional gravity 
theory. It satisfies the 5D Ricci-flat field equation Rab — 0. In this paper, two branes are embedded into the bulk. 
One brane is at y — where the standard matter lives. The other brane is at y = y%, where y\ is the thickness of 
bulk. It is the basal topological structure of this black string space. The usual 4D effective cosmological constant A 
is considered to be induced from the 5D Ricci-flat space. One should notice that A is considered as a free parameter. 
The distance between black hole horizon and cosmological horizon is shorten with bigger A. In this metric it has 
nothing to do with the value referred from current cosmological observation. So, if the value of A increases, the Nariai 
black hole is inevitably arisen in its last fate. 

2. For the well known Nariai case, we have solved the scalar field around it in the 5D SdS black string space. The 
fifth dimensional component L{y) is presented. Taking into account the classical field theory, one know that standard 
model fields (such as fermions, gauge bosons, Higgs fields) are confined on a (3 + 1) dimensional hypersurface (3-brane) 
without accessing along the transverse dimensions. In order to stabilize two branes, the scalar field is led in. According 
to "standing wave" condition, the fifth dimensional equation can be solved. The scalar field gets its maximum on the 
brane y = and get an extremeness value on the other brane y — y\. So the spectrum of parameter f2 is broken 
into two parties, one is quantum f2 n and the other is continuous one. The quantum spectrum is illustrated in Fig. [1] 
It is clear that the extreme Nariai black hole decays more acutely than usual case [241 ]. Furthermore, the quantum 
phenomenon emerges distinctly in the waves. One can see those according to the effective potential (|24l) . Fig. [3] and 
Fig. [3 

3. Because of the singular f(r) — in the metric, potential (flU]) vanishes both near black hole horizon and 
cosmological horizon. Then the Schrodinger wavelike equation (J23J) reduces to a solvable one (p?5|) . Obviously, according 
to the real scalar field, we get the boundary conditions ^Ff\ . Eq. (|2"3")) describes one dimensional transmission of waves 
scattering by potential barrier. In order to solve this equation, we adopt a useful tangent approximation Q to unite 
radial coordinate r and tortoise coordinate x. So Eq. (|2"3")l , effective potential (j2"4"]) and boundary condition l[27p 
constitute a full boundary value problem. Because of the complicated potential (I24[) and fitting function ([25)) , we only 
give numerical solution. By used the replacement, the continuous potential is switched into a square barrier. With 
the classical method of the square barrier, the reflection and transmission coefficients (R, T) are obtained naturally. 
The result is presented briefly in Tabic fl] 
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